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ABSTRACT 

We present the results of time-resolved spectroscopy of 13 0-type stars in 
the Gas OB6 stellar association. We conducted a survey for radial velocity vari¬ 
ability in search of binary systems, which are expected to be plentiful in young 
OB associations. Here we report the discovery of two new single-lined binaries, 
and we present new orbital elements for three double-lined binaries (including 
one in the multiple star system HD 17505). One of the double-lined systems is 
the eclipsing binary system DN Gas, and we present a preliminary light curve 
analysis that yields the system inclination, masses, and radii. We compare the 
spectra of the single stars and the individual components of the binary stars with 
model synthetic spectra to estimate the stellar effective temperatures, gravities, 
and projected rotational velocities. We also make fits of the spectral energy dis¬ 
tributions to derive E{B — V), R = Av/E{B — V), and angular diameter. A 
distance of 1.9 kpc yields radii that are consistent with evolutionary models. We 
hnd that 7 of 14 systems with spectroscopic data are probable binaries, consis¬ 
tent with the high binary frequency found for other massive stars in clusters and 
associations. 
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Subject headings: binaries: close — stars: early-type — stars: individual (DN Cas, 
BD-|-60°497, HD 17505, HD 17520, BD-|-60°594) — open clusters and associa¬ 
tions: individual (Cas OB6, IC 1805, IC 1848) 


1. Introduction 

The frequency of spectroscopic binaries (SBs) in star forming regions has been a topic 
of interest for many years. Mason et ah (1998) review the fraction of binary systems found 
among O stars in clusters/associations, the held, and runaway systems. The results from 
their large survey agree with previous studies (e.g. Garmany, Conti, & Massey 1980; Levato et 
ah 1991) and indicate that a large fraction of O stars in clusters and associations are binaries 
while fewer held stars show multiplicity and very few runaway stars have companions. The 
values from Mason et al. (1998) are broken down into spectroscopic and resolved binaries. For 
O stars in clusters and associations they give a total of 75% as binaries; 61% of the binaries 
are spectroscopic and 41% are resolved (with some obvious overlap). The spectroscopic 
binary group includes, however, a signihcant number of cases labeled as SBl? or SB2? 
(systems suspected of radial velocity variability or showing double lines), so the fraction of 
conhrmed spectroscopic binaries is then 34% among O stars in clusters and associations. 

We have conducted a radial velocity survey of 0-type stars in the Cas OB6 stellar 
association in search of spectroscopic binary systems. The Cas OB6 region is one of the 
largest complexes in the Perseus arm of the Galaxy and is located at a distance of about 
2.3 kpc (Garmany & Stencel 1992; Massey, Johnson, & DeGioia-Eastwood 1995; Foster & 
Routledge 2003). The warm dust and ionized gas in this region extends over 150 pc along 
the arm and is dominated by the W3/W4/W5 chain of H II regions and the HB3 supernova 
remnant (Carpenter, Heyer, & Snell 2000; Terebey et al. 2003). The open cluster IC 1805 
(OCL 352 = Melotte 15) is situated at the western end in the W4 region and contains 8 
0-type stars (Massey et al. 1995). The three brightest of these, HD 15558, 15570, and 15629, 
are probably the ionization source that powers the blowout of gas in the Galactic chimney 
(Normandeau, Taylor, & Dewdney 1996; Terebey et al. 2003) and the emission arc above 
W4 (Reynolds, Sterling, & Haffner 2001). The open cluster IC 1848 (OCL 364) is found 
in the eastern part of the association in the W5 region. Its primary source of ionizing flux 
is the multiple 0-star system HD 17505 (Stickland & Lloyd 2001). In addition to these 
two clusters. Carpenter et al. (2000) hud evidence of 19 other young clusters that are still 
embedded in their primordial clouds. 
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Our knowledge of the binary star population among the young, massive stars of Gas OB6 
is quite limited at present. There are a number of radial velocity investigations of the sample 
but the overall small number of measurements has made the detection of binary stars difficult 
(Hayford 1932; Underhill 1967; Ishida 1970; Conti, Leep, & Lorre 1977; Liu, Janes, & Bania 
1989, 1991; Rauw & De Becker 2004). Spectroscopic orbits are available for only three 
systems, HD 15558 (Garmany & Massey 1981), HD 16429 (McSwain 2003), and BD+60°497 
(Rauw & De Becker 2004), while a fourth star, BD+60°470 = DN Gas, is a known eclipsing 
binary (Frazier & Hall 1974). However, results for 0-stars in other clusters and associations 
suggest that many more binaries may exist (Mason et ah 1998; Garcia & Mermilliod 2001). 
It is important to determine the extent of the binary star population in Cas OB6 in order 
to estimate accurately the ionizing flux of the key 0-type stars and to assess the role of 
gravitational encounters with binary stars that might lead to ejection of runaway stars (Gies 
& Bolton 1986; Hoogerwerf, de Bruijne, & de Zeeuw 2001). 

Here we present a reconnaissance of the radial velocity variations observed in 13 bright 
stars in Cas OB6 that we obtained over a 5 night run with the KPNO 4 m Mayall telescope 
in 2003 October. We also present the hrst radial velocity study of the eclipsing binary sys¬ 
tem DN Cas, which we combine with published photometry to determine the component 
masses. We discuss our observations and reduction techniques in §2. The radial velocity 
measurements and the derivation of orbital elements are summarized in §3. We use a to¬ 
mography algorithm to separate the spectra of the binary star components, and we describe 
a comparison of the stellar spectra with model spectra in §4 that we use to estimate the 
basic stellar parameters. We examine in §5 the spectral energy distributions to determine 
reddening, ratio of total-to-selective extinction, and angular size for each target, and we 
discuss the probable distance to Cas OB6. Details about the individual objects follow in §6. 
Finally we summarize in §7 our derived fraction of multiple and single systems for Cas OB6 
in the context of past results. 


2. Spectroscopic Observations and Rednctions 

The blue spectra were obtained with the RC spectrograph at the Mayall 4 m telescope 
at Kitt Peak National Observatory from 2003 September 30 through 2003 October 4. The 
spectra were made with the BL 380 grating (1200 grooves mm“^, blazed at 4500 A in second 
order) using a BG39 order sorting hlter. The detector was the T2KB CCD, a 2048 x 2048 
device with 24 pm square pixels. The spectra recorded wavelengths between 4180 and 
4940 A, although internal vignetting occurs for wavelengths red-ward of about 4800 A. The 
reciprocal dispersion of the spectra is approximately 0.37 A pixeU^, and the resolving power 
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is R = A/AA = 5700. Exposure times were usually a few minutes in length, and each stellar 
spectrum was preceded and followed by a He-Ne-Ar lamp exposure for careful wavelength 
calibration. Most of the resulting spectra have a signal-to-noise ratio in excess of 100 pixel“^. 
We usually obtained at least two spectra of each target on each of hve consecutive nights. 
We also made numerous bias and flat held frames each night. A list of the names of targets 
appears in Table 1. 

We also obtained an additional 32 spectra of HD 17505A with the KPNO 0.9 m Coude 
Feed telescope in runs in 2000 October and 2000 December (Boyajian et al. 2005). The 
spectra were made with the long collimator, grating B (in second order with order sorting 
hlter OG550), camera 5, and a Ford 3072 x 1024 CCD (F3KB) detector. This arrangement 
produced a resolving power, R = A/AA = 9530, and covered a range of 829 A around Ha. 
Exposure times were 10 - 15 minutes, and the resulting spectra have a S/N = 200 in the 
continuum. We also observed the rapidly rotating A-type star, ( Aql, which we used for 
removal of atmospheric water vapor lines. Each set of observations was accompanied by 
numerous bias, hat held, and Th Ar comparison lamp calibration frames. 

The spectra were reduced using standard routines in IRAF^ to create continuum rectihed 
versions of each spectrum. Then all the spectra of a given target were transformed to a 
common heliocentric frame in log A form and collected in a spectrum matrix of wavelength 
and time for subsequent analysis. The rectihcation of the vignetted portion of the spectrum 
was generally problem-free, so we used the entire spectral range in the following analysis. 


3. Radial Velocities and Orbital Elements 

We found that the spectra of 10 of the 13 targets were composed of one set of spectral 
lines, and for these targets we obtained radial velocities by hnding the extremum point of the 
cross-correlation of each spectrum with a spectrum template (omitting those spectral regions 
containing interstellar features). We initially obtained velocities using the hrst spectrum 
as the template, and then we shifted each spectrum to the reference frame of the first 
spectrum and co-added all the spectra to form a high S/N template. The hnal relative 
velocities were then found by cross-correlation with this high S/N template spectrum. These 
relative velocities were then transformed to absolute velocity by adding the velocity measured 
in the template spectrum by fitting parabolae to the lower portions of selected spectral 


"‘IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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lines. Rest wavelengths were taken from the list of Bolton & Rogers (1978). The hnal 
velocities for these 10 stars are listed in Table 2 (which appears only in the electronic version), 
and the mean velocity < V). > for each star is given in column 3 of Table 1. Note that 
the adopted absolute velocity of the template spectrum depends critically on what subset 
of lines is measured because of the presence of signihcant line-to-line velocity differences. 
These differences mainly result from line formation at differing heights in an expanding 
stellar atmosphere (Hutchings 1976; Bohannan & Garmany 1978; Gies & Bolton 1986). We 
generally selected high excitation lines that form deep in the atmosphere where the outflow 
velocity is smallest (see discussion of the individual targets in §6). We list in column 4 of 
Table 1 — the standard deviation of the measured velocity among different lines in the 

template spectrum. 

We used a version of the statistical F-test to determine if a star was a velocity variable 
object. Since we usually obtained at least two spectra every night of each object and since we 
do not expect large variations in velocity over the time span of a few hours, we estimated the 
average measurement error from the mean of the absolute value of the differences within a 
night divided by \/2. This spectrum-to-spectrum error estimate (t(s — s) is given in column 

5 of Table 1. We then compared this to the overall standard deviation of the velocities 
a (column 6 of Table 1), and set the criterion for variability as a/a{s — s) > 3.3, which 
corresponds to the <1% probability that both variances are drawn from the same population 
(for 9 degrees of freedom = 10 observations minus 1 mean value). We found that three of 
the ten single-lined targets were velocity variable according to this criterion (see the hnal 
column of Table 1). The stars HD 17520 and BD-|-60°594 are probably long period binary 
systems while the very luminous star HD 15570 exhibits variations related to atmospheric 
structure huctuations (see the discussion in §6). 

The spectra of the targets DN Gas and BD-|-60°497 both showed evidence of line- 
doubling in some observations, and the spectra of HD 17505A displayed three components 
in the best separated cases. We were able to use a tomography algorithm to isolate the sta¬ 
tionary line component in HD 17505A (§6.5), and after removal of the stationary component 
the spectra were treated in the same way as for the other double-lined systems. We measured 
velocities for these double-lined binaries in a two-step process (Gies et al. 2002). Preliminary 
radial velocities were found using the multiple Gaussian htting technique within the IRAF 
SPLOT routine. Then preliminary orbital elements were determined that provided starting 
estimates of the velocities for both components in each spectrum. We then formed template 
spectra for both components by interpolating in the 0-star synthetic spectrum grid of Lanz 

6 Hubeny (2003) for a temperature and gravity appropriate for each star, and these were 
convolved with broadening functions to account for rotational and instrumental broadening 
(§4). Then we obtained velocities for specihc spectral lines by a non-linear least squares £t 
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of the observed profile with the two template spectra prohles co-added with a pre-dehned 
monochromatic flux ratio (details are given for each case in § 6 ). We took care to form an 
average velocity for each spectrum based on the same set of selected lines in order to avoid 
spurious variations introduced by line-to-line velocity differences. The resulting velocities 
are given below in §6 (Tables 4, 6 , and 8 ). 

We determined orbital elements for these three double-lined binaries using the non¬ 
linear, least-squares htting routine of Morbey & Brosterhus (1974). The starting value for 
the orbital period was obtained with the period search technique of Scargle (1982) as modihed 
by Horne & Baliunas (1986) (or from photometry in the case of DN Cas). We first solved for 
the best £t period and then we formed separate elliptical solutions for both the primary and 
secondary components. The derived eccentricity in the case of DN Cas was not significantly 
different from zero (Lucy & Sweeney 1971), but in the other two targets the eccentricity 
was significant and the same within errors in the solutions for the primary and secondary 
components. We formed average values of the eccentricity e, longitude of periastron uj, and 
time of periastron T that were then fixed in solutions to obtain the semiamplitude K and 
systemic velocity 7 for the primary and secondary components individually (since differing 
atmospheric expansion rates may lead to different systemic velocities). The final orbital 
elements are given in §6 (Tables 5, 7, and 9). 


4. Stellar Parameters 

We co-added our spectra to form high S/N versions for use in spectral classification, 
temperature and gravity estimation, and measurement of the projected rotational velocity 
V sin i. We used a simple shift-and-add algorithm to form an average spectrum in the case 
of the single-lined targets. However, we relied upon the Doppler tomography algorithm 
described by Bagnuolo et ah (1994) to reconstruct the individual component spectra of the 
multiple component systems. The tomography algorithm is an iterative scheme that uses 
the derived orbital velocity shifts and monochromatic flux ratio to estimate the spectra of 
each component. We ran the algorithm for 50 iterations with a gain of 0.9, although the 
results are insensitive to both parameters. The final, continuum rectified spectra of all our 
targets (except HD 17505; see §6.5) and all their components are shown in Figures 1 and 2 
in order of decreasing effective temperature. 

We estimated the spectral classification of each target (Table 1) by comparing their 
spectra with standard star spectra in the atlas of Walborn & Fitzpatrick (1990). The cor¬ 
responding stellar effective temperature Tgs was determined by comparing the equivalent 
widths of various He II/He I line ratios with those measured in the synthetic spectral library 
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for 0-stars by Lanz & Hubeny (2003). In most cases, we simply assumed a surface gravity of 
loggf = 4.0 as a suitable value for main sequence stars, but this approach was inadequate in 
several cases of evolved stars where the H Balmer line wings are narrower and indicative of a 
lower gravity. We experimented with a number of model synthetic spectra to best estimate 
the gravity in these cases. This was accomplished through comparing models by eye with the 
data in order to hud an appropriate match for surface gravity. The most difficult spectrum 
to match was that of the very luminous star HD 15570, which displays clear evidence of 
the atmospheric expansion into a stellar wind that is not included in the static atmosphere 
formulation adopted by Lanz & Hubeny (2003). We estimate that our htting errors are 
approximately ±1000 K in temperature and ±0.3 dex in logg, and the results are listed in 
Table 1. 

We estimated the projected rotational velocity V sini by comparing the observed He I 
and He II line widths with those found by convolving the model prohles from Lanz & Hubeny 
(2003) with instrumental and rotational broadening functions. The instrumental broadening 
of our spectra was represented by a Gaussian with FWHM = 52 km s“^. The rotational 
broadening function (Gray 1992) was computed assuming a spherical star with a linear limb- 
darkening coefficient determined from Wade & Rucinski (1985). A grid of theoretical prohles 
was computed for a range of V sini, and a minimization was used to determine the best 
ht to the data. The results and errors are listed in Table 1. Note that this approach does 
not account for variations in the intensity prohle across the disk, nor the non-spherical shape 
and gravity darkening found in rapid rotators (Townsend, Owocki, & Howarth 2004). These 
shortcomings are most important for the two rapid rotators in our sample, BD±60°513 and 
BD±60°594, and we have probably underestimated the true rotational broadening in these 
cases. 


5. Radii from Fits of the Spectral Energy Distributions 

All of our targets are considered to be members of the Gas OB6 association based 
upon their position in the sky and placement in the Hertzsprung-Russell diagram (HRD) 
(Voroshilov et al. 1985; Garmany & Stencel 1992; Walborn 2002). The proper motion study 
of IG 1805 made by Vasilevskis, Sanders, & van Altena (1965) conhrms membership for hve 
targets but curiously indicates a somewhat discrepant proper motion for HD 15558. We will 
assume that all the targets share a common distance, so their observed fluxes are directly 
related to their relative stellar radii. Since we have determined temperature and multiplicity 
for each target, we have all parameters needed to estimate their radii as a function of distance 
by comparing their observed and model flux distributions, and we can then check the distance 



by comparing the derived and model radii (as a function of temperature). The angular 
diameter of the limb darkened disk dm (in units of radians) is found by the inverse-square 
law: 


/a (observed) _ ^ _ 1 2 

Fa (emitted) * 4 ^ 


( 1 ) 


where the ratio of the observed and emitted fluxes (reduced by the effects of interstellar 
extinction Aa) depends on the square of the ratio of stellar radius F* to distance d. All our 
targets have a large reddening, so we need to £t the relation above over the entire observed 
spectrum in order to account reliably for the wavelength dependent shape of the extinction 
curve. We adopted the extinction curve law from Fitzpatrick (1999) that is a function of the 
reddening E{B — V) and the ratio of total-to-selective extinction R = Ay/E{B — V). We 
then fit the observed flux distribution of each target to determine E{B — V), R, and angular 
diameter dm- 


We compiled data on the observed fluxes from the UV to the near-IR. There are low 
dispersion International Ultraviolet Explorer (lUE) spectra available (SWP and LWR cam¬ 
eras) for all of the targets except DN Cas, HD 17520, HD 237019, and BD-|-60°586, and 
we formed binned versions of the UV spectral fluxes from these. We found Johnson UBV 
photometry for all the targets (plus Cousins RI in a few cases), and these magnitudes were 
transformed to fluxes using the calibration of Colina, Bohlin, & Castelli (1996). We also 
included Stromgren photometry for 7 of the 13 targets using the calibration of Gray (1998). 
All of the stars are included in the 2MASS All-Sky Catalog of Point Sources (Cutri et al. 
2003), and we converted the JHK magnitudes to fluxes using the calibration of Cohen, 
Wheaton, & Megeath (2003). 

The model fluxes were selected from the grid of Lanz & Hubeny (2003) based upon 
the values of Teff and logg in Table 1. These model fluxes are based upon plane-parallel, 
non-LTE, line blanketed atmospheres that do not include the effects of stellar wind mass 
loss. However, the inclusion of winds generally alters the spectral energy distributions only 
at wavelengths much larger and much smaller than those considered here (Martins, Schaerer, 
& Hillier 2005), so wind effects are insignificant in this context. We co-added model fluxes 
using the observed flux ratios for the spectroscopic binaries and using the observed magnitude 
difference for close visual binaries (details are given in next section for individual stars). Note 
that in the two single-lined binary cases, HD 17520 and BD-|-60°594, we simply assumed 
that the unseen companion contributes 5% of the total flux (a compromise between the 
expectation that the mass ratio is not too extreme and that brighter companions would have 
visible spectral features in our observations). We have no spectral information for some of 
the companions, and in these cases we assumed that the companion was a main sequence 
star with a temperature corresponding to the observed flux ratio (from the calibration in 
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Table 4 of Martins et al. 2005). We set a lower limit of = 27500 K and loggf = 4.0 
for the faint companions, which corresponds to the cool boundary in the grid of Lanz & 
Hubeny (2003). This simplihcation is acceptable since any detected companion is probably 
hot enough that the optical and near-IR regions correspond to the common Rayleigh-Jeans 
portion of the flux distribution. 

We compared the observed and model flux distributions and made a grid search £t to 
find the values of E{B — V) and R that produced the lowest residuals (by weighting 
the observations approximately equally between the UV, optical, and near-IR bands). The 
normalization term in the hnal £t yields the angular diameter. Our results are gathered in 
Table 3 which lists the star name, home cluster, angular diameter (in micro-arcsec), derived 
radius (see below), E{B — V), and R. The acccompaning error estimates for E{B — V) 
and R are derived from the point in grid search where xl attains a value higher by 1 than 
normalized minimum value, and the errors in diameter are based upon the scatter in flux 
distribution fit and the errors propagated from uncertainty in Tefj, E{B — V), and R. The 
model hts were generally excellent, but they were less satisfactory for the spectra of the 
three hottest stars, where it was difficult to find a consistent match of the far- and near- 
UV fluxes. We also found the JHK fluxes were too large relative to the optical fluxes for 
HD 17520, which we suspect is due to an IR excess associated with the Be disk emission of 
the companion (Porter & Rivinius 2003), so the fit was based on the optical fluxes alone in 
this case. 

The stellar radius is found by 


R* = 107.43 Rq Old d (2) 

where the distance d is given in kpc. The distance estimates for IC 1805 and IC 1848 range 
from 1.9 kpc (Ishida 1970) to 2.4 kpc (Guetter & Vrba 1989) (with a scatter of ±10% among 
recent estimates). We selected a distance of d = 1.886 kpc from the compendium of open 
cluster data by Loktin et al. (2001) because this low end value appears to give the best match 
between the observed and theoretical stellar radii and provides good agreement with radii of 
similar stars in eclipsing binaries (Gies 2003). The resulting radii are listed in Table 3 and 
are plotted against stellar temperature in Figure 3. There is a clear and consistent trend 
for the hotter main sequence stars to be larger, and the observations match the predictions 
for main sequence stars from the atmosphere models of Martins et al. (2005) (shown as a 
dotted line in Fig. 3). We also show in Figure 3 the (Teg, R*) loci for stars with ages of 1.00 
and 3.16 Myr from the evolutionary models of Lejeune & Schaerer (2001), which correspond 
to the age range for IG 1805 given by Massey et al. (1995). We see that most of the stars 
fall within the expected region in this diagram for this age range. However, the positions 
of several stars indicate that there is a real dispersion in the ages among this sample. For 
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example, BD+60°501 appears to have a relatively small radius for its temperature, indicating 
it is the youngest object in our sample. The derived stellar radii are directly proportional 
to the assumed distance, and, for example, if we adopted the larger distance d = 2.34 kpc 
advocated by Massey et ah (1995) then the stars would be approximately 24% larger, a 
size that would be difficult to reconcile with the evolutionary models of Lejeune & Schaerer 
(2001). We show below that our derived radii for DN Cas agree well with those determined 
independently from an analysis of the eclipsing light curve. 


6. Results for Individual Objects 
6.1. Stars With Constant Radial Velocity 

BD+60°501.- The two Balmer lines Hy and H/5 have the lowest measured velocities in 
the template spectrum, so they were omitted from the average (i.e., based only on the He I 
and He II lines). Our observed mean velocity of —57.9 ± 1.1 km s“^ is reasonably consistent 
with the mean from the work of Rauw & De Becker (2004) of —49.9 ± 2.5 km s“^ who also 
hnd the star to be a constant velocity object. 

HD 15558.- There is evidence here for a velocity progression (lower velocity among low 
excitation lines), so we used only the lines He II AA4199,4541 to set the radial velocity of the 
template spectrum. There are no obvious variations over the 5 night run, but the star may 
be a long period, single-lined binary (Garmany & Massey 1981). According to the recent 
orbital elements from Stickland & Lloyd (2001), the orbital period is 439.64±0.26 d and the 
predicted velocity at the mean time of our observations is —21.8 km s“^ (near the peak of 
the radial velocity curve), which is 27.7 km s“^ higher than our observed mean of —49.5±2.9 
km s“^. This discrepancy could be caused by differences in measurement techniques or errors 
in the preliminary orbital solution of Stickland & Lloyd (2001). A long-term study by M. 
De Becker and G. Rauw (2005, private communication) has recently conhrmed the binary 
status of this star and will provide an improved orbital solution. Therefore, we list this star 
as “constant (SBl)” in Table 1 denoting both a lack of short term line variability in our 
observations as well as the recent conhrmation of binarity. There is a visual companion at a 
separation of 9'.'9 that is 2.6 mag fainter (Perryman 1997). 

HD 15570.- There is a clear velocity progression between spectral lines with the low 
excitation lines (H, He I) appearing at a much lower velocity than the high excitation lines 
(He II, N V). There is also a correlation with line width in the sense that stronger, broader 
lines have lower velocity. We adopted a radial velocity for the template spectrum from the 
lines He II A4199 and N V AA4604,4620. The night-to-night scatter between velocities is 
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small except on the last night when the measured velocity is ~ 10 km s“^ more positive 
than the mean from the other nights. Curiously, the Of emission lines do not show any shift 
in velocity on the hnal night. We suspect that the change was related to an atmospheric 
fluctuation rather than orbital motion. Underhill & Gilroy (1990) also consider this star to 
be a velocity variable. De Becker, Rauw, & Eenens (2005) hnd no evidence of orbital motion 
and discuss rotational modulation as a possible source of line variability. 

HD 15629.- The low excitation Balmer and He I lines have systematically low velocities, 
so the template mean was adopted from the lines He II AA4199,4541,4686. We find no 
evidence of velocity variability. Underhill & Gilroy (1990), on the other hand, flag this star 
as possibly variable, but given the presence of line-to-line velocity differences, some of the 
velocity differences found by various authors are probably due to the line sample adopted 
(and not to orbital motion). De Becker et ah (2005) also hnd no evidence of orbital motion 
for this star. 

BD+60°513.- This star has a broad-lined spectrum, and we selected only the strongest 
lines to measure the template mean velocity (Hy, H/3, He I AA4471,4921, He II AA4541,4686). 
We conhrm the lack of short term velocity variability reported by Rauw & De Becker (2004), 
and our mean velocity of —58.6 ± 2.8 km s“^ is probably consistent with their mean value, 
-44.3 ± 10.5 km s'h 

BD-|-62°424.- The H/5 line had the lowest velocity measured, so we omitted the H 
Balmer lines and weak lines in forming the template mean velocity. The star appears to be 
radial velocity constant over the run, and our mean velocity, —34.3 ± 1.7 km s“^, is probably 
consistent with that from Abt et al. (1972) of —43.5 ± 6.3 km s“h 

HD 237019.- All the strong lines were used to form the template mean velocity. The 
star has a constant radial velocity over the run. Our mean velocity of —43.4 ±1.1 km s“^ 
is higher than the only other value in the literature of —68 km s“^ from Seyfert & Popper 
(1941). 

BD±60°586.- We omitted H/5 and He I A4921 from the template mean calculation, since 
both lines had signihcantly lower velocities. The star was radial velocity constant during the 
run. Our mean velocity, —48.1 ± 1.0 km s“^, agrees with that found by Abt et al. (1972), 
—42.3 ± 7.1 km s“\ but is lower than the value reported by Gonti et al. (1977), —27.9 ± 1.2 
km s“^. The star is a member of a visual multiple system (ADS 2194) with companions at 
separations of 772 (AU = 1.7) and 1370 (AU = 4.7). 
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6.2. Single-lined Spectroscopic Binaries 

HD 17520.- The H/9 line has a blue-shifted emission component which also appears 
weakly in Hy. The Balmer lines and weak lines were omitted in the calculation of the 
template mean velocity. The velocities show a progressive increase over the duration of 
the run, reaching values similar to those observed by Conti et ah (1977) and Liu et al. 
(1989). Walter (1992) showed how Ho emission developed between 1985 and 1991 leading 
to a spectral classihcation as a Be-type star (probably its status in 2003 as well since the 
emission is so apparent in H/5). This star is also the central object in a visual multiple 
system (ADS 2165) with a close companion at a separation of 0'.'335 and with a magnitude 
difference of 0.52 mag (Perryman 1997; Walborn 2002). Our spectra show evidence of only 
one set of lines (for the 09 V component), and unfortunately, we do not know if the Be 
emission originates in the 0-star or its nearby companion (which could have broad shallow 
lines as observed in many Be stars). We note that over the course of the run the H/5 emission 
appeared to be almost stationary while the absorption lines shifted red-ward by ~ 17 km s“^. 
The simplest interpretation is that the 0-star component is a long-period, single-lined binary, 
while the Be emission originates in the fainter component of the close visual pair. There is 
another star at a separation of 1078 that is 3.2 mag fainter in K (Cutri et al. 2003). 

BD-|-60°594.- We used the strong lines to set the template mean velocity for the broad- 
lined spectrum of this star. The first measured velocity is close to the value obtained by 
Conti et al. (1977), —50.5 ±5.6 km s“\ but the subsequent velocities decline monotonically 
by ~ 49 km s“^ over the next 4 days. This suggests that the star is a single-lined binary 
with a period of ~ 20 d and a semiamplitude of ~ 50 km s“^. 


6.3. The Double-Lined Eclipsing System DN Cas 

DN Cas has been known as an eclipsing binary for some time and was classified as 
an Algol type binary with a period of 1.155479 d by Hoffmeister (1947). Hiltner (1956) 
classihed the spectrum as 08 V. A later photometric study by Frazier & Hall (1974) used 
UBV photoelectric photometry to rehne the orbital period to 2.310955 d, twice that of the 
original value. Unfortunately, the comparison star used by Frazier & Hall (1974) turned 
out to be itself an eclipsing binary, compromising their hnal photometry. Davidge (1980) 
made another photometric study, and based upon the curvature near the maxima of the 
light curve he reclassihed DN Cas as a /3 Lyrae type system. Interestingly no spectroscopic 
observations of DN Cas have appeared in the literature, though Davidge (1980) claimed 
that such observations were underway. Additionally, the only photometry that has been 
published, other than in phase-folded hgures, is the compromised data of Frazier & Hall 
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(1974). This overlooked opportunity to obtain masses for this massive binary is rectified here, 
where we present the first full orbital solution based on our spectroscopy and photometry 
from the literature. 

We determined velocities for both components using non-linear least-sqnares fits of two 
templates with temperatnres and rotational velocities given in Table 1 and with a monochro¬ 
matic flnx ratio in the bine of F 2 /F 1 = 0.41 ± 0.10. Onr hnal radial velocities are given in 
Table 4 and the spectroscopic orbital elements are given in Table 5. The period was adopted 
from the eclipse timing observations collected by Kreiner, Kim, & Nha (2001). Here we adopt 
the time of minimnm light (primary snperior conjnnction) to set the zero-point of orbital 
phase and epoch Tmin, and our derived epoch agrees with the current photometric epoch of 
minimnm light from Kreiner et al. (2001) of HJD 2,452,501.928 ±0.007. Eccentric hts offered 
no signihcant improvement over circular hts (Lucy & Sweeney 1971), so circular elements 
are given in Table 5. The radial velocities and orbital solntions are shown in Fignre 4. 

The only pnblished photometric data of DN Cas are those from Frazier & Hall (1974) 
who made UBV photoelectric observations between 1971 October and 1972 December. Un- 
fortnnately, their comparison star, HD 14817, is itself an eclipsing binary as mentioned 
above. However, the light curve of HD 14817 from Hipparcos (Perryman 1997) is hat out¬ 
side of eclipse, so we can retain the resnlts from Frazier & Hall (1974) that are based on 
outside-of-eclipse times for HD 14817. These measurements are plotted against orbital phase 
in Fignre 5, and because there remain signihcant phase gaps in this restricted set, onr light 
curve htting results should be considered preliminary. 

We calculated a model light curve using the Wilson-Devinney code (Wilson & Devinney 
1971; Wilson 1990). Most of the model parameters are set by the spectroscopic orbital 
solntion (Table 5), the adopted temperatnres and gravities from the spectra (Table 1), the 
associated limb darkening parameters (Van Hamme 1993), and the monochromatic hux ratio 
(which sets the ratio of the mean stellar radii). In addition, we set the bolometric albedo 
and gravity darkening exponents to 1.0, which are valid approximations for stars with fnlly 
radiative envelopes. We then ht the observed V-band observations by varying the inclination 
i and the stellar radii. The ratio of stellar radii was set using the observed monochromatic 
flux ratio and model snrface flux ratio (set according to the temperatures and gravities in 
Table 1) 

R2/R1 = [(/A(observed)2//A(observed)i) x (FA(emitted)i/FA(emitted)2)]°'^ = 0.72 ±0.09. 

If the stars are rotating synchrononsly, then the ratio of projected rotational velocities also 
yields the ratio of radii as 


R2/R1 = iy sini) 2 / iy sini)i = 0.72 ± 0.15 
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and the good agreement between these two estimates snggests that synchronous rotation 
has been attained in this close binary. We found that the radii obtained from a £t of the 
spectral energy distribution and distance of 1.9 kpc (Table 3) led to an excellent fit of the 
light curve for an orbital inclination of i = 79°2 ± 1?1 (Fig. 5), and the derived stellar masses 
are Mi = 19.2 ± 0.8 Mq for the 08 V primary and M 2 = 13.9 ± 0.5 Mq for the BO.2 V 
secondary (Table 5). The mass of the primary is consistent with its spectral type for the 
calibration given by Martins et ah (2005), and the mass of the secondary agrees well with 
the mass for its spectral type according to the relations given by Hanson, Howarth, & Conti 
(1997). Both stars are well within their respective Roche lobes. 

There is a companion to DN Cas in the 2MASS catalogue at a separation of 573 and at 
a position angle of 89°. However, inspection of the 2MASS images at this position suggests 
that this object is probably not a real star since its point spread function looks very sharp 
and unlike those of nearby stars. Furthermore, we see no evidence of such a tertiary in 
our observations. Since the near east-west alignment of the putative object and DN Cas 
is the same as the slit orientation for our KPNO 4 m spectra, the tertiary should appear 
in our spectra with a separation from DN Cas of 7.7 pixels in the spatial dimension. The 
average 2MASS magnitude difference is 2.13 mag (for FT, 77), so that the flux ratio of tertiary 
to binary is 0.14 in the near-IR. However, our blue spectra show no sign of a companion 
brighter than Am = 5, so if this companion is real, it is probably a very red foreground 
star that can be ignored in our analysis. We found that the introduction of a third light 
contribution of 14% in the R-band light curve results in an increase in the inclination by 3°, 
but since the actual flux contribution of a tertiary is probably much smaller than this, any 
changes in our results due to the presence of a tertiary are probably comparable to or less 
than the errors quoted in Table 5. 


6.4. The Double-Lined System BD-|-60°497 

The spectra of BD-|-60°497 show clear evidence that this star is a double-lined spec¬ 
troscopic binary. In their recent paper, Rauw & De Becker (2004) give the first orbital 
solution for BD-|-60°497, and they find an orbital period of 3.96 days. We measured ve¬ 
locities using the scheme outlined in §3 by matching template synthetic profiles for models 
with temperatures and gravities given in Table 1, together with an estimated flux ratio 
of F 2 /-F 1 = 0.35 ± 0.08. We note, however, that in some of the lines (for example. He I 
AA4387,4471) we found that the approaching, blue-shifted component (primary and sec¬ 
ondary) often appeared to be weaker than expected based upon this fixed flux ratio (Rauw 
& De Becker 2004 also observed a similar variation). We doubt that this variation has a 
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significant effect on our measurements, which appear in Table 6. 

We combined our velocities with those of Rauw & De Becker (2004) in making the 
orbital solution. A number of alias periods remain viable, but we obtained the best £t with 
a period close to that obtained by Rauw & De Becker (2004) (see Table 7). However, unlike 
Rauw & De Becker (2004), we found that a significant eccentricity is present, and we made 
full elliptical fits for both components. Otherwise our results are in reasonable agreement 
with their solution (compared in Table 7). The velocities and orbital fit are shown folded on 
the derived ephemeris in Figure 6. Assuming masses typical of their spectral classes (Martins 
et al. 2005), the system probably has an inclination of i ~ 47°, so we do not expect eclipses 
to occur. 


6.5. The Multiple Star System HD 17505 

The star system HD 17505 dominates the center of the cluster IC 1848, and there are 
seven visual companions noted by Mason et al. (1998) at distances ranging from 2'.'1 to 12470. 
Multiplicity in the central object HD 17505A has been discussed by Conti & Alschuler (1971), 
Walborn (1973), Howarth et al. (1997), and Walborn (2002) and it is known to contain at 
least three stars since the single International Ultraviolet Explorer spectrum shows evidence 
of three Doppler shifted components (Stickland & Lloyd 2001). Walborn (2002) pointed out 
that the multiplicity of the star is the reason it appears so bright in the HRD. We obtained 
10 spectra of HD 17505A over the five night duration of the KPNO 4 m run, and we found 
that the lines were strongly blended in all but a few spectra corresponding to the maximum 
Doppler separation of the binary motion. However, we also obtained an additional 32 red 
spectra of the target in 2000 with the KPNO 0.9 m Coude Feed telescope (§3), and these 
turned out to be particularly helpful in finding the orbital elements. 

Initial inspection of these higher dispersion red spectra showed that all three compo¬ 
nents were cleanly resolved in the He I profiles observed near orbital quadrature phases and 
indicated an orbital period of approximately 8 days. We measured radial velocities by first 
making unconstrained three-Gaussian fits of the least blended feature. He I A7065, in the 
quadrature phase spectra. We obtained a preliminary orbital solution from this subset of 
spectra for the two moving components, and then we used this solution to estimate the 
Doppler shifts for all the Coude Feed observations. Then we made three-component Gaus¬ 
sian fits of He I A7065 in all the red spectra using these starting estimates for the Doppler 
shifts and constraining the widths to be the mean values found from the fits of the quadra¬ 
ture phase spectra. A new orbital solution was obtained, and the solution was used to make 
a three component reconstruction of the individual component spectra using the Doppler 
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tomography algorithm (Bagnuolo et ah 1994). We subtracted the reconstructed spectrum of 
the stationary component from each spectrum and renormalized them to end up with a set 
corresponding to the spectra of the double-lined binary alone. We then measured the radial 
velocities using the non-linear, least-squares method by matching template spectra to the 
prohles of Ha, He 1 A6678/He 11 A6683, and He 1 A7065. The template spectra for the two 
components were taken from their respective tomographic reconstructions {F 2 /F 1 = l.OiO.l) 
rather than using the model prohles from Lanz & Hubeny (2003), because we were concerned 
that some of these red features are not well represented by the synthetic models. This ap¬ 
proach worked well except in one conjunction phase observation (from HJD 2,451,892.743) 
where the lines were too badly blended for reliable measurement. Our hnal radial velocities 
from the red spectra for the close binary are collected in Table 8. 

The estimate for the orbital period from the red spectra alone was P = 8.572 ± 0.003 d, 
and the error was sufficiently small that we could determine the cycle number and approx¬ 
imate orbital phase corresponding to the times of the KPNO 4 m observations in 2003 
October. We measured the line width of the He 1 A4471 feature in these blue spectra to hnd 
that maximum velocity separation occurred at HJD 2,452,915.469 ±0.12, and we compared 
this to the time of the same orbital phase for the 2000 Coude Feed observations to arrive 
at a revised period estimate of P = 8.5710 ± 0.0008 d. The hnal orbital solution is given 
in Table 9 and the radial velocity curve appears in Figure 7. The predictions from this ht 
for the time of the lUE observation are ±128 and —184 km s“^, in reasonable agreement 
with the observed velocities of ±142 and —165 km s“^ (Stickland & Lloyd 2001), but we 
did not include the lUE measurements in the hnal solutions because of possible systematic 
diherences between velocities measured in the UV and optical. 

We made tomographic reconstructions of both the blue and red spectra using the or¬ 
bital solution in Table 9 and assuming hux contributions of 30%, 30%, and 40% for Aal, 
Aa2, and Ab, respectively (based upon the relative equivalent widths of the He 1 A6678 ± 
He 11 A6683 feature). Both the blue and red spectral reconstructions are plotted in Figure 
8 . We determined radial velocities of the stationary component Ab of —30 ± 3 and —38 ± 5 
from the blue and red reconstructed spectra, respectively. Once again, the optical velocities 
are systematically more negative than the UV determined value for Ab of —13.5 km s“^ 
(Stickland & Lloyd 2001). Our results are summarized in a schematic representation of the 
system hierarchy in Figure 9. Note that the Coude Feed spectra were made with a 275 wide 
slit in an approximately north-south orientation, so the B component (at a separation of 
2715 and a position angle of 93° with a magnitude difference of AHp = 1.66 mag; Perryman 
1997) may be slightly mixed in with the light from the Ab component for the red recon¬ 
structed spectrum. However, the 4 m spectra were made with an east-west slit orientation, 
so it was possible to separate the A and B components spatially along the slit. A single 
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extracted spectrum of the B component is shown in Figure 8 along with the tomographic 
reconstructions of the Aal, Aa2, and Ab components. The spectrum of the B component 
is clearly different and indicates a cooler effective temperature than the stars of component 
A, and, thus, the spectroscopic triple we detect is clearly related to A alone. The hottest 
and brightest of all the components is Ab, and its spectral type agrees with that obtained 
by Walborn (1972) for the entire system. 

The stars of the Aal-Aa2 central binary are far enough apart that no eclipses should be 
seen, and indeed none are observed (Perryman 1997). For example, if we adopt the masses 
and radii associated with 0-type main sequence stars from Martins et ah (2005), then the 
expected system inclination should be i = 65° (based upon the Msin^i values in Table 9) 
while eclipses would be observed only if i > 73° (based upon the semimajor axis in Table 9). 
The orbital period of B around the triple A can be estimated as ~ 27, 000 y if we take the 
projected angular separation as the semimajor axis and the distance is 1.9 kpc. The orbital 
period of the Aa-Ab central triple is probably much less. The system was unresolved in the 
speckle survey of Mason et ah (1998), and if we take their resolution limit of 35 mas as the 
projected semimajor axis then the orbital period should be < 61 y. The target is clearly 
worthy of continued spectroscopic and interferometric observation to search for evidence of 
such orbital motion. 


7. Summary 

We have monitored thirteen O stars in the Cas OB6 association in search of multiple star 
systems. Five of the thirteen were found to have radial velocity variations related to orbital 
motion. Of these five, we have determined orbital solutions for three double-lined systems: 
the eclipsing binary DN Cas, BD-|-60°497, and the spectroscopic triple HD 17505A. The 
remaining two binaries, BD-|-60°594 and HD 17520, are single-lined binaries with periods 
longer than 5 days. In addition to these hve, one of our targets, HD 15558, is probably a long 
period binary (Garmany & Massey 1981; Stickland & Lloyd 2001). One other bright 0-star 
in Cas OB6 is HD 16429, which was fonnd by McSwain (2003) to be a spectroscopic triple. 
If we include all of these objects in the known sample of observed 0-stars in the Cas OB6 
association, then we hnd that 7 of the 14, or 50%, are spectroscopic binaries. This value is 
well within the range fonnd by Mason et al. (1998) for spectroscopic binaries among stars in 
clusters and associations. 

We note that two systems, HD 17505 and DN Cas, are of special interest and would 
beneht from follow-up observations. We are planning to obtain new three-color photometry 
of the eclipsing binary DN Cas in order to make a better £t of the light curve and to determine 
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an accurate mass of the 08 V primary star. The multiple system HD 17505 is a fascinating 
system with at least four O stars (with a total mass close to 100 Mq) which are apparently 
gravitationally bound. High resolution studies via interferometry to separate the Aa and Ab 
components would produce a more accurate kinematic model for the system and would help 
us better understand stellar interactions in these systems and their role in the formation and 
ejection of massive stars. 
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Table 1. Target Object Summary 


Object 

Name 

Spectral 

Class. 

< Vr 
(km s~^) 

ail — V) 
(km s“^) 

(j(s — s) 
(km s“^) 

CF 

(km s~^) 

Tff 

(kK) 

log g 
(cm s“^) 

V sin 2 
(km s“^) 

Binary 

Status 

DN Cas A 

08 V 

-52 


4.9 

149.3 

34.2 

4.0 

158 ± 20 

SB2 

DN Cas B 

B0.2 V 

-38 


5.3 

207.3 

29.7 

4.0 

113 ± 18 

SB2 

BD+60°497 A 

06 V 

-43 


2.2 

126.8 

37.5 

4.0 

183 ±9 

SB2 

BD+60°497 B 

08 V 

-81 


3.6 

166.7 

33.0 

4.0 

171 ± 19 

SB2 

BD+60°501 

06.5 V 

-58 

9.4 

0.5 

1.1 

36.8 

4.0 

173 ± 10 

constant 

HD 15558 

05 HI(f) 

-49 

5.3 

1.8 

2.9 

39.5 

3.7 

195 ± 11 

constant (SBl) 

HD 15570 

04 If+ 

-55 

0.9 

0.8 

4.2 

42.9 

3.5 

123 ± 10 

atmospheric var. 

HD 15629 

05 V((f)) 

-60 

8.1 

1.2 

2.2 

41.5 

4.0 

106 ±6 

constant 

BD+60°513 

07.5 Vn 

-59 

6.7 

1.5 

2.8 

34.8 

4.0 

259 ±5 

constant 

BD+62°424 

06.5 V((f)) 

-34 

3.7 

0.8 

1.7 

36.2 

3.8 

103 ± 12 

constant 

HD 17505 Aal 

07.5 V((f)) 

-9 


8.6 

115.3 

36.3 

4.0 

90± 10 

SB2 

HD 17505 Aa2 

07.5 V((f)) 

-44 


7.2 

118.2 

36.5 

4.0 

120 ± 10 

SB2 

HD 17505 Ab 

06.5 HI((f)) 

-38 

8.2 



37.5 

3.8 

80± 15 

constant 

HD 17520 

09 V 

-54 

2.7 

0.7 

6.0 

33.2 

4.0 

80± 12 

SBl + Be 

HD 237019 

08 V 

-43 

8.1 

0.7 

1.1 

34.0 

4.0 

128 ±5 

constant 

BD+60°586 

07.5 V 

-48 

4.7 

1.0 

1.0 

35.4 

4.0 

95 ±5 

constant 

BD+60°594 

08.5 Vn 

-86 

8.3 

1.0 

17.7 

33.5 

4.0 

285 ± 16 

SBl 


Systemic velocities for the binary stars are given in Tables 5, 7, and 9. 
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Table 2. Radial Velocity Measurements for Constant Velocity and SBl Stars 


Star 

Name 

Date 

(HJD-2,450,000) 

(km s“^) 

BD+60°501 

2912.891 

-56.6 

BD+60°501 

2912.968 

-56.4 

BD+60°501 

2913.888 

-57.2 

BD+60°501 

2913.965 

-57.9 

BD+60°501 

2914.890 

-57.4 


Note. — Table 2 is published in its entirety in 
the electronic edition of the Astrophysical Jour¬ 
nal. A portion is shown here for guidance re¬ 
garding its form and content. 







Table 3. Spectral Energy Distribution Fits 


Object 

Name 

Gluster 

Membership 

(^LD 

(10“® arcsec) 

p a 

(Re) 

E{B - V) 
(mag) 

R 

{Av/E{B - H)) 

DN Gas A .... 

IG 1805 

36(1) 

7.4 (7) 

0.97 (1) 

3.13 (1) 

DN Gas B .... 

IG 1805 

26 (3) 

5.3 (8) 

0.97 (1) 

3.13 (1) 

BD+60°497 A 

IG 1805 

49 (4) 

10.0 (13) 

0.89 (2) 

2.94 (6) 

BD+60°497 B 

IG 1805 

32 (5) 

6.5 (11) 

0.89 (2) 

2.94 (6) 

BD+60°501 .. 

IG 1805 

36 (3) 

7.2 (9) 

0.76 (2) 

3.12 (5) 

HD 15558 A .. 

IG 1805 

81 (10) 

16.4 (26) 

0.84 (3) 

3.07 (9) 

HD 15558 B .. 

IG 1805 

28 (4) 

5.6 (9) 

0.84 (3) 

3.07 (9) 

HD 15570 .... 

IG 1805 

95 (16) 

19.3 (37) 

1.01 (2) 

3.13 (7) 

HD 15629 .... 

IG 1805 

57 (7) 

11.6 (19) 

0.76 (3) 

3.09 (12) 

BD+60°513 .. 

IG 1805 

43 (4) 

8.7 (12) 

0.83 (3) 

3.00 (5) 

BD+62°424 .. 

IG 1805 

49 (3) 

9.9 (11) 

0.75 (2) 

3.02 (4) 

HD 17505 Ab . 

IG 1848 

56 (6) 

11.4 (16) 

0.72 (3) 

2.83 (9) 

HD 17505 Aal 

IG 1848 

49 (8) 

10.0 (18) 

0.72 (3) 

2.83 (9) 

HD 17505 Aa2 

IG 1848 

49 (8) 

10.0 (18) 

0.72 (3) 

2.83 (9) 

HD 17505 B .. 

IG 1848 

44 (6) 

8.9 (15) 

0.72 (3) 

2.83 (9) 

HD 17520 Aa . 

IG 1848 

38 (11) 

7.7 (24) 

0.57 (5) 

3.10 (120) 

HD 17520 Ab . 

IG 1848 

14 (4) 

2.8 (9) 

0.57 (5) 

3.10 (120) 

HD 17520 B .. 

IG 1848 

34 (10) 

7.0 (22) 

0.57 (5) 

3.10 (120) 

HD 17520 G .. 

IG 1848 

14 (4) 

2.8 (9) 

0.57 (5) 

3.10 (120) 

HD 237019 ... 

IG 1848 

37(1) 

7.4 (7) 

0.77 (1) 

3.20 (6) 

BD+60°586 A 

IG 1848 

45(3) 

9.1 (11) 

0.58 (3) 

3.27 (32) 

BD+60°586 B 

IG 1848 

25 (2) 

5.1 (7) 

0.58 (3) 

3.27 (32) 

BD+60°594 A 

IG 1848 

35 (4) 

7.2 (11) 

0.68 (2) 

2.88 (10) 

BD+60°594 B 

IG 1848 

10 (2) 

2.0 (5) 

0.68 (2) 

2.88 (10) 


'^Radius derived for a distance of 1.886 kpc. 
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Table 4. Radial Velocity Measurements for DN Cas 


Date 

(HJD-2,450,000) 

Orbital 

Phase 

Vi 

(km s~^) 

(O - C)i 
(km s“^) 

V2 

(km s“^) 

(O - C)2 

(km s“^) 

2912.849 

0.814 

143.9 

-1.4 

-318.3 

-7.1 

2912.932 

0.850 

118.0 

-3.8 

-282.5 

-3.8 

2913.003 

0.881 

99.4 

4.6 

-238.3 

2.9 

2913.825 

0.237 

-255.4 

2.7 

244.1 

-4.2 

2913.834 

0.241 

-257.6 

0.8 

245.4 

-3.4 

2913.839 

0.243 

-259.7 

-1.1 

247.4 

-1.6 

2913.844 

0.245 

-256.1 

2.6 

246.6 

-2.5 

2913.927 

0.281 

-255.9 

-1.0 

241.1 

-2.7 

2914.004 

0.314 

-241.9 

0.1 

221.2 

-4.7 

2914.814 

0.665 

133.1 

0.4 

-297.8 

-4.0 

2914.943 

0.721 

154.9 

-3.7 

-337.5 

-7.7 

2915.025 

0.756 

154.1 

-8.0 

-352.3 

-17.8 

2915.656 

0.029 

-109.8 

-23.1 

4.7 

-5.9 

2915.883 

0.127 

-200.8 

-1.5 

161.9 

-4.8 

2915.931 

0.148 

-221.7 

-4.6 

182.5 

-9.0 

2916.015 

0.184 

-239.8 

1.4 

229.9 

5.1 

2916.747 

0.501 

-59.8 

-13.1 

-47.3 

-2.4 

2916.756 

0.505 

-31.6 

10.1 

-82.0 

-30.1 

2916.792 

0.521 

-17.3 

3.4 

-66.2 

14.7 

2916.797 

0.523 

-14.6 

3.6 

-56.4 

28.2 

2916.805 

0.526 

-10.9 

3.0 

-55.6 

34.9 

2916.809 

0.528 

-18.8 

-7.5 

-59.1 

35.1 

2916.818 

0.532 

3.7 

10.3 

-102.2 

-1.6 

2916.822 

0.534 

8.7 

12.6 

-98.6 

5.7 

2916.846 

0.544 

17.7 

8.4 

-123.8 

-1.1 

2916.851 

0.546 

14.4 

2.5 

-128.6 

-2.3 

2916.943 

0.586 

64.2 

4.3 

-195.9 

-3.2 

2917.018 

0.618 

92.2 

-1.9 

-247.1 

-6.8 






Table 5. Orbital Elements for DN Gas 


Element 

Value 

P (days) . 

2.310950 (2) 

T^in (HJD-2,450,000) 

2915.579 (2) 

e . 

0 

^ (°) . 


Ki (km s“^) . 

210.5 (20) 

K 2 (km s“^) . 

292.0 (38) 

M 2 /M 1 . 

0.721 (12) 

7 (km s“^) . 

-48.3 (14) 

7 (km s“^) . 

-42.7 (27) 

in . 

79.2 (11) 

Ml (Me) . 

19.2 ( 8 ) 

M 2 (Me) . 

13.9 (5) 

Ri {R&) . 

7.4 ( 6 ) 

R 2 {Rq) . 

5.4 ( 8 ) 

i?i/i?i(lobe) . 

0.78 ( 6 ) 

R2/R2i}ohe) . 

0.66 ( 10 ) 

® {Rq) . 

23.3 (2) 

r.m.s.i (km s“^) . 

7.4 

r.m.s .2 (km s“^) . 

14.1 


Note. — Numbers in parentheses 
give the error in the last digit quoted. 
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Table 6. Radial Velocity Measurements for BD+60°497 


Date 

(HJD-2,450,000) 

Orbital 

Phase 

Vi 

(km s“^) 

(O - C)! 
(km s“^) 

^2 

(km s“^) 

(O - C)2 

(km s“^) 

2912.871 

0.151 

-214.8 

0.6 

136.3 

-6.1 

2912.953 

0.171 

-224.3 

-6.5 

145.4 

-0.1 

2913.874 

0.404 

-100.0 

2.7 

30.7 

34.7 

2913.946 

0.422 

-66.7 

22.4 

-18.9 

2.7 

2914.823 

0.644 

70.7 

8.0 

-216.0 

2.7 

2914.826 

0.644 

65.0 

2.0 

-217.6 

1.7 

2914.832 

0.646 

63.6 

-0.3 

-217.5 

2.9 

2914.836 

0.647 

60.1 

-4.3 

-230.5 

-9.5 

2914.844 

0.649 

64.5 

-0.9 

-220.8 

1.4 

2914.848 

0.650 

68.6 

2.6 

-222.9 

0.2 

2914.857 

0.652 

69.8 

2.6 

-219.4 

5.1 

2914.862 

0.654 

63.9 

-3.8 

-228.9 

-3.6 

2914.965 

0.680 

74.6 

-5.4 

-238.0 

3.1 

2915.825 

0.897 

41.3 

-3.9 

-194.0 

1.9 

2915.829 

0.898 

39.8 

-4.4 

-197.5 

-2.8 

2915.833 

0.899 

36.6 

-6.6 

-198.4 

-5.0 

2915.842 

0.901 

38.9 

-2.0 

-191.8 

-1.3 

2915.846 

0.902 

38.0 

-2.0 

-194.2 

-5.0 

2915.952 

0.929 

16.9 

6.6 

-169.0 

-18.3 

2916.831 

0.151 

-210.9 

4.6 

142.0 

-0.5 

2916.835 

0.152 

-206.3 

9.3 

148.4 

5.7 

2916.859 

0.158 

-214.1 

2.5 

139.0 

-5.0 

2916.863 

0.159 

-216.2 

0.6 

145.6 

1.4 

2916.932 

0.176 

-212.9 

5.0 

139.0 

-6.7 

2917.011 

0.196 

-219.1 

-2.8 

132.5 

-11.1 
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Table 7. Orbital Elements for BD+60°497 


Element Rauw & De Becker (2004) This paper 


P (days) . 3.96 (9) 3.95863 (21) 

T (HJD-2,450,000) ■ ■ ■ 2916.233 (15) 

To (HJD-2,450,000) 2935.98 (15) 

e . 0 0.156 (19) 

cu (°) . ■■■ 100 ( 11 ) 

Ki (km s"^) . 172 (12) 159.9 (22) 

K 2 (km s-i) . 221 (15) 207.7 (29) 

M 2 /M 1 . 0.78 ( 8 ) 0.770 (15) 

7 1 (km s-i) . -54 (9) -53.8 (17) 

72 (km s-i) . -69 (12) -67.5 (23) 

Ml sin^i (Mq) . 13.9 (25) 11.1 (5) 

M 2 sin^i (Mq) . 10.9 (20) 8.6 (4) 

a sin i (T©) . 30.7 (19) 28.39 (28) 

r.m.s.i (km s“^) ... ■ ■ ■ 10.9 

r.m.s .2 (km s“^) ... ■ ■ ■ 14.4 


Note. — Numbers in parentheses give the error in the last digit 
quoted. 
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Table 8. Radial Velocity Measurements for HD 17505Aa 


HJD 

(-2,450,000) 

Orbital 

Phase 

Vi 

(km s“^) 

(O - C)i 
(km s“^) 

V2 

(km s“^) 

(O - C)2 

(km s~^) 

1817.887 

0.772 

-181.1 

5.3 

143.7 

5.2 

1818.890 

0.889 

-184.6 

-2.8 

136.3 

2.6 

1819.873 

0.004 

-77.9 

-1.3 

18.9 

-6.9 

1820.887 

0.122 

63.4 

-0.2 

-124.7 

-6.7 

1821.864 

0.236 

148.2 

15.8 

-190.4 

-1.9 

1822.865 

0.353 

120.9 

3.7 

-171.9 

1.1 

1823.814 

0.464 

49.2 

0.0 

-105.9 

-2.9 

1823.924 

0.476 

35.7 

-3.4 

-95.7 

-2.8 

1824.825 

0.582 

-45.5 

4.8 

1.0 

2.0 

1824.950 

0.596 

-54.8 

8.4 

16.6 

4.5 

1830.842 

0.284 

142.9 

7.5 

-187.7 

3.8 

1830.946 

0.296 

136.5 

2.5 

-180.0 

10.1 

1888.797 

0.045 

-29.1 

-4.5 

-34.0 

-6.5 

1889.746 

0.156 

91.2 

-2.1 

-148.5 

-0.0 

1892.743 

0.506 





1893.774 

0.626 

-77.2 

11.8 

25.5 

-13.0 

1894.771 

0.742 

-168.5 

4.5 

131.3 

6.6 

1894.845 

0.751 

-183.3 

-6.0 

136.3 

7.2 

1895.721 

0.853 

-196.3 

-1.8 

144.7 

-2.1 

1895.811 

0.864 

-184.3 

7.5 

144.7 

0.6 

1896.651 

0.962 

-132.7 

-8.0 

76.7 

1.4 

1896.784 

0.977 

-103.8 

4.1 

45.3 

-12.6 

1897.652 

0.078 

6.0 

-10.1 

-71.1 

-1.8 

1897.784 

0.094 

21.0 

-12.9 

-87.4 

0.1 

1898.656 

0.196 

114.5 

-3.8 

-177.2 

-3.2 

1898.793 

0.212 

126.0 

0.8 

-173.9 

7.2 

1899.658 

0.312 

133.2 

2.5 

-182.4 

4.3 

1899.791 

0.328 

124.9 

-1.5 

-182.2 

0.2 

1900.653 

0.429 

69.7 

-4.9 

-136.9 

-7.7 

1900.784 

0.444 

60.6 

-3.2 

-124.9 

-6.7 

1901.637 

0.543 

-32.7 

-15.8 

-14.1 

21.2 

1901.770 

0.559 

-27.6 

3.0 

-24.4 

-3.0 
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Table 9. Orbital Elements for HD 17505Aa 


Element 

Value 

P (days) . 

8.5710 (8) 

T (HJD-2,450,000) 

1862.696 (16) 

e . 

0.095 (11) 

^ (°) . 

252 (6) 

Ki (km s“^) . 

166.5 (18) 

K 2 (km s“^) . 

170.8 (18) 

M 2 /M 1 . 

0.975 (15) 

7 i (km s“^) . 

-25.8 (12) 

72 (km s“^) . 

-26.3 (12) 

Ml sin^i {Mq) . 

17.1 (6) 

M 2 sin^i (ATq). 

16.6 (6) 

a sin i (i?©) . 

56.8 (4) 

r.m.s.i (km s“^) .. 

7.0 

r.m.s .2 (km s“^) .. 

7.0 


Note. — Numbers in parentheses 
give the error in the last digit quoted. 
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4200 4400 4600 4800 

Wavelength (A) 


Fig. 1.— The observed mean spectra sorted by temperature. Here we show the seven hottest 
stars with the hottest at the top. The spectra have been continuum normalized and shifted 
vertically for clarity. Note that interstellar features (such as the 4428 A diffuse interstellar 
band) are cut out of the tomographic reconstructions but are present in the other spectra. 
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4200 4400 4600 4800 


Wavelength (A) 


Fig. 2.— The observed spectra sorted by temperature. Here we show the seven coolest stars 
with the coolest at the bottom. The spectra have been continuum normalized and shifted 
vertically for clarity. 
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Ten (kK) 


Fig. 3.— The stellar radii (for an assumed distance of 1.9 kpc) plotted versus effective 
temperature {diamonds = V, squares = I, III luminosity classes). The solid and dashed lines 
show the predicted radii for ages of 1.00 and 3.16 Myr, respectively (Lejeune & Schaerer 
2001). The dotted line indicates the relationship given for main sequence stars by Martins 
et ah (2005) (see their Table 4). 







1 


plotted 
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Orbital Phase (2452916.233(15) + 3.95863(21)) 


Fig. 6.— Our radial velocity data {circles) with the data from Rauw & De Becker (2004) 
{triangles) for BD+60°497. The data is folded on the given ephemeris and plotted with the 
calculated model {line). The solid points denote radial velocities of the primary and the 
open points represent secondary star velocities. 
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Diagram showing the hierarchy of the HD 17505 system. 



